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Filmwise  condensation  of  steam  on  externally-finned  tubes 
is  a  very  complex  process.  Recent  experiments  have  shown 
that  enhancement  ratios  (ratio  of  steam-side  heat-transfer 
coefficient  to  that  of  a  smooth  tube  having  the  same  diameter) 
exceeded  the  area  enhancement  produced  by  the  fins.  Moreover, 
the  enhancement  ratios  for  fully  flooded  tubes  exceeded  the 
values  predicted  by  a  simple,  one-dimensional  conduction 
model  by  a  factor  of  2  to  4 .  A  new  two-dimensional  conduction 
model  was  developed,  which  showed  that  the  one-dimensional 
model  overpredicted  the  two-dimensional  results  for  high 
conductivity  tube-metals  such  as  copper  by  as  much  as  13%. 

The  two-dimensional  model  also  showed  that  variations  in  fin 
thickness  or  spacing  can  result  in  an  overprediction  by  the 
one-dimensional  model  of  the  two-dimensional  results  by  as 
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I .  INTRODUCTION 


A.  BACKGROUND 

-^Because  the  size  and  weight  of  main  propulsion  and 
auxiliary  systems  are  inversely  proportional  to  the  ease  in 
which  heat  or  energy  is  exchanged/  a  major  thrust  of  the 
research  and  development  program  of  the  U.S.  Navy  is  toward 
the  design  and  development  of  equipment  that  perform  at 
higher  efficiency  with  a  reduction  in  size  and  weight.  In 
particular,  one  area  of  great  interest  is  the  reduction  in 
size  and  weight  of  steam  condensers^ 

Current  equipment  is  of  the  ^re  or  smooth  tube  type  of 
heat  exchanger  in  which  the  ste^  is  condensed  on  the  outside 
of  the  tubes  (the  shell  side) .  '^he  heat-transfer  effective¬ 
ness  is  governed  by  the  amounf.  of  surface  and  the  overall 
resistance,  to  the  flow  of  heat^  This  overall  resistance,  in 
turn,  is  the  sum  of  the  vapor-side,  tube  wall,  and  water-side 
thermal  resistances.  Reduction  of  any  or  all  of  these 
resistances  results  in  an  increase  in  heat-transfer  performance 
with  an  associated  reduction  in  condenser  surface  area  which 
translates  into  decreased  size. 

Both  the  vapor-side  and  tube-side  thermal  resistances  can 
be  considered  as  the  reciprocal  of  the  so-called  hS  product 
where  h  is  the  coefficient  of  heat  transfer  and  S  is  the 


surface  area  through  which  the  heat  flows.  Because  of  the 
curved  surface  of  the  exchanger  tubing,  the  vapor- "space" 


and  tube-side  surfaces  are  not  equal  and  an  adjustment  must 
be  made  to  the  tube-side  heat-transfer  coefficient  to  allow 
for  a  design  on  the  basis  of  the  vapor  or  shell-side  surface. 

It  is  a  fact  that  the  vapor-side  hs  product  (reciprocal 
of  resistance)  can  be  as  high  as  two  or  three  times  the  tube- 
side  hS  product.  Thus,  it  would  appear  reasonable,  in  an 
attempt  to  decrease  the  overall  resistance,  to  attack  and  try 
to  increase  the  hS  product  on  the  tiabe  side.  Unfortunately, 
little  can  be  done  in  the  arena  of  increasing  h  significantly 
on  either  side  although  the  vapor  side  has  greater  possibili¬ 
ties  in  this  regard  and,  by  its  very  location,  not  much  can 
be  done  with  regard  to  the  increase  of  the  tube-side  surface. 

It  is  for  the  foregoing  reasons  that  recent  efforts  have 
focused  on  the  use  of  externally-finned  tubes  in  steam  con¬ 
densers  to  increase  the  flow  of  heat  by  increasing  the  vapor- 
side  hS  product  [Refs.  1,2,3]. 

For  the  past  few  decades,  it  has  been  accepted  that  the 
use  of  externally-finned  tubes  would  be  inappropriate  mainly 
because  of  the  liquid  retention  (flooding)  that  occurs  between 
fins  (Figure  1.1).  This  flooding  is  a  result  of  the  relatively 
high  surface  tension  of  water  (at  50°C  water  has  a  surface 
tension  of  0.068  N/m,  while  refrigerants  possess  values  about 
0.015  N/m).  One  of  the  unavoidable  results  of  flooding  is  the 
considerably  lower  heat  transfer  through  the  flooded  portion 
than  that  through  the  unflooded  portion.  For  the  case  of  an 
externally  finned  copper  tube  with  steam  as  the  condensate, 


the  copper  has  a  thermal  conductivity  600  times  greater  than 
that  of  water  (385  W/m-K  for  copper  compared  with  0.6  W/m*K 
for  water) .  Due  to  its  overall  importance  in  the  heat- 
transfer  scheme,  the  flooding  of  the  lower  portion  of  the  tube 
by  condensate  has  been  studied  theoretically  and  experimentally 
by  several  investigators  [Refs.  4, 5, 6, 7, 8].  These  studies  have 
shown  that  the  flooding  angle  depends  mainly  on  the  fin 
spacing  and  the  ratio  of  surface  tension  to  density  of 
condensate  as  given  by  the  following  equation  for  the  case  of 
rectangular-section  fins  [Ref.  5] : 


Recent  work  by  Yau  et  al.  [Ref.  9]  and  Wanniarachchi  et  al.  [Ref.  1] 
has  shown  that  even  in  the  cases  with  considerable  flooding, 
the  heat-transfer  performance  exceeds  that  expected  on  the 
basis  of  area  increase  due  to  the  fins. 

In  1948,  Beatty  and  Katz  [Ref.  10]  developed  the  first  model 
to  predict  the  condensation  heat-transfer  coefficient  for 
horizontal,  externally-finned  tubes.  Their  model  ignores 
the  existence  of  flooding.  Although  their  model  agrees  to 
within  10  percent  for  condensing  fluids  possessing  low  surface 
tensions,  it  tends  to  be  rather  inaccurate  for  high-surface- 
tension  fluids,  such  as  water.  Several  other  investigators, 

Rudy  and  Webb  [Ref.  2]  and  Owen  et  al.  [Ref.  8]  attempted  refinement 
of  the  Beatty  and  Katz  model  but  still  were  unable  to  predict 
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the  heat-transfer  coefficient  to  a  reasonable  degree  of 
accuracy. 

Rudy  and  Webb  used  the  Beatty  and  Katz  model  for  the 
unflooded  portion,  while  neglecting  heat  transfer  through  the 
flooded  portion.  With  this  modification,  Rudy  and  Webb  under¬ 
predicted  data  for  low-surface-tension  fluids  by  up  to 
50  percent.  Owen  et  al.  attempted  to  improve  the  Rudy  and 
Webb  model  by  allowing  heat  transfer  through  the  flooded 
portion.  In  order  to  compute  heat  transfer  through  the 
flooded  portion,  they  assumed  parallel  conduction  paths 
through  the  fins  and  the  condensate  that  is  flooded  in  the 
interfin  space  as  shown  in  Figure  1.2.  With  this  modification, 
they  improved  the  predictions  slightly;  underprediction  for 
low-surface-tension  fluids  decreased  to  about  40  percent. 
However,  as  shown  by  Honda  and  Nozu  [Ref.  11], both  of  these 
models  overpredicted  the  steeim  data  of  Yau  et  al.  by  a  factor 
of  up  to  3.  The  fact  is  that  the  use  of  Beatty  and  Katz 
model,  which  does  not  take  surface  tension  forces  into  consi¬ 
deration  for  the  unflooded  portion,  may  be  responsible  for 
this  large  discrepancy. 

The  most  recent  and  accurate  attempt  at  modelling  the 
problem  was  accomplished  by  Honda  and  Nozu  [Ref.  11].  Their 
resulting  solution,  a  fourth-order  differential  equation  for 
the  film  thickness,  was  solved  nvimerically  subject  to  reasonable 
boundary  conditions.  An  agreement  to  within  ±20%  was  achieved 
by  their  model  in  predicting  data  experimentally  obtained 
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Fiaure  1.2.  One-Dimensional  Parallel  Conduction  Kodel 


for  11  fluids  and  22  tubes.  However,  the  disagreement  with 
steam  condensation  data  was  as  high  as  40%. 

To  understand  the  condensing  heat  transfer  performance 
more  clearly,  it  is  convenient  to  define  the  enhancement  ratio 
as  proposed  by  Wanniarachchi  et  al.  [Refs. 1,2],  The  enhancement 
ratio  is  defined  as  the  ratio  of  the  steam-side  heat-transfer 
coefficient  of  a  finned  tube  to  that  of  a  smooth  tube  with  the 
same  diameter  as  the  finned  tube  root  diameter  at  the  same 
heat  flux. 

The  present  investigator  has  been  mainly  motivated  by  the 
heat  transfer  through  a  finned  tube  that  is  completely  flooded 
by  condensate.  As  shown  by  Winniarachchi  et  al.  [Ref .  2],  a  finned 
tube  with  a  fin  spacing  of  0.5  mm,  and  a  height  and  thickness 
of  1.0  mm  and  1.0  mm  respectively  showed  a  considerable 
improvement  in  heat-transfer  performance  even  though  this  tube 
was  completely  flooded  by  condensate.  In  fact,  this  tube 
showed  an  enhancement  ratio  of  about  2.5  under  vacuum  (85  mmHg) 
and  4.5  at  atmospheric  pressure.  These  ratios  even  exceeded 
the  area  enhancement  (A^/A^),  which  was  2.11.  If  one¬ 
dimensional  conduction  is  assumed  (see  Owen  [Ref.  8] ) ,  the  expected 
enhancement  ratio  was  computed  to  be  about  0.8.  This  unex¬ 
pectedly  high  enhancement  ratio  may  be,  at  least  in  part, 
explained  by  two-dimensional  conduction  and  convection  in 
the  interfin  space. 


OBJECTIVES 


The  main  objectives  of  this  thesis  are: 

Develop  a  two-dimensional  numerical  model  of  a  horizontal 
finned  tube  fully  flooded  with  steam  condensate.  The 
program  using  the  model  shall  be  interactive  to  allow 
different  cases  of  fin  thickness,  fin  spacing,  and 
tube-metal  thermal  conductivity  to  be  studied. 

Using  a  thermal  analyzer  computer  program,  determine 
the  temperature  distribution  plots  and  heat  flux  for 
each  case. 

Compute  enhancement  ratios  for  each  one  of  the  cases 
studied  and  compare  these  ratios  with  those  computed 
by  using  a  one-dim.ensional  model, 

Com.pare  the  one  and  two-dimensional  enhancement  ratios 
and  comment  on  the  importance  of  the  role  of  two- 
dimensional  heat  transfer  in  the  overall  heat-transfer 
performance , 

Determine  the  value  and  usefulness  of  the  two-dimensional 
model  as  an  aid  in  the  prediction  of  enhancement  ratios 
for  various  cases. 


Consider  the  three-dimensional  Laplace  steady-state 


equation  for  the  spatial  temperature  distribution  in  a 
configuration 


0 


(2.1) 


It  may  be  solved  analytically,  subject  to  certain  boundary 
conditions,  and  the  solution  is  called  an  analytical  solution 
The  finite-difference  solution  is  arrived  at  numerically  and 
applies  to  a  network  of  points  called  nodes,  with  the 
temperatures  at  the  nodes  being  estimated.  The  numerical 
solutions  obtained  by  analytical  and  finite-difference 
methods  can  be  made  to  agree  with  any  desired  accuracy. 

Consider  a  grid  of  nodes.  Applications  of  equation 
(2.1)  at  a  particular  node  whose  temperature  is  T^  allov/s 
estimation  of  the  second  derivatives  in  terms  of  the 
temperature  T^  and  the  six  related  temperatures  at  the 
adjacent  nodes  of  the  grid.  In  Figure  2.1  an  arbitrary 
function  T  =  f (x)  is  shown  which  considers  only  the  trace 
of  T  =  f2(3<)  even  though  equation  (2.1)  shows  that  there 
will  also  be  a  functional  relationship  between  T  and  y, 
namely,  T  =  f^ (y)  and  between  t  and  z,  namely,  T  = 


22 


If  y  and  z  are  held  constant,  the  first  derivative 
of  equation  (2.1)  may  be  approximated  by  the  change  in  T 
divided  by  the  change  in  x: 


AT  +  AT 

AXj^  +  AX2 


(T2  -  T^)  +  (Tq  -  T^)  T2  - 

Ax^  +  AX2  ~  Ax^  + 


(2.2) 


An  estimate  of  the  second  derivative  is  obtained  by 
estimating  the  first  derivative  in  each  of  the  two  intervals 
and  dividing  the  difference  of  the  first  derivatives  by  the 


total  interval  over  which  the  first  derivative  estimates 
apply.  The  total  interval  will  be  the  distance  between  the 
centers  of  the  two  intervals,  Ax^^  and  AX2.  Hence 


AT2/AX2  -  ATj^/AXj^ 
2  AX2  +  2 


If  Ax^  =  AX2  =  AX, 


(T2-To)/Ax  -  (T^-Tj^)/Ax 


T,  +  T_  -  2T^ 
12  O 


A  similar  development  using  T  =  ^2^^^  v;ith  temperatures 
T,  and  T.  vields 


and  for  T  =  £3(2)  with  temperatures  and 


9^T 

82^ 


Tj.  +  T-  -  2T 
5 _ 6 _ o 

Az 


(2.5) 


so  that  equation  (2.1)  may  be  approximated  by 


T,  +  T-  -  2T  T,  +  T.  -  2T  T-  +  T,  -  2T 

1  2  ^ - 2^_5 - ^ °  =  0(2.6) 


Ax 


Ay 


Az 


and  if  Ax  =  Ay  =  Az,  the  result  is 


T,  +  T_  +  T,  +  T.  +  T-  +  T,  -  6T  =  0 

1  2  3  4  5  6  o 


Further  justification  of  the  result  indicated  by 
equation  (2.3)  comes  from  a  consideration  of  a  Taylor  Series 
Expansion.  Recall  that  the  Taylor  Series  Expansion  is  given 
by 


f(x)  =  ^(x-a)O  .  ^(x-a)  .  ^^(x-a)2 


f (a) 


f"(a) 


f^(x-a)-"  +  ...  +  ---;^(x-a) 

•  1 1  • 


n 


(2.7) 


Two  cases  will  be  examined.  First  for  f (x)  =  T(x) 


at  x  =  x^,  then 


T(x=x^)  =  T^ 


T(x=x^)  =  T^ 


and  with  these  in  equation  (2.7) 


o  dx 


1 

d^T  2 

.  1 

d^T 

Ax  +  2 

— ^  A  X 

+  6 

x_ 

dx  X 

dx 

+  ^  AX^'  +  .  .  . 

dx"'  x„ 
o 


(2.8) 


Then  for  f (x)  =  T(x)  at  x  =  X2 


T(x=X2)  =  T2 
T(x=x^)  = 
a  =  X 
x-a  =  -AX 


With  these  in  equation  (2.7) 


dT 

^2  =  ^o  -  d3^ 


1  d^Tl 
Ax  +  TT  - o 


1  d'*!  ,  4 

+  24  77 

dx  X 

o 


'  dx'  x„ 


(2.9) 


VThen  equations  (2.8)  and  (2.9)  are  added,  one  obtains 


T-  +  T, 


=  „  .  5^(Ax)'  .  1  ^(4x)«  + 


(2.10) 


If  Ax  is  small,  then  all  terms  above  the  second  order 


in  Eq.  (2.9)  may  be  considered  negligible.  Thus  upon 
rearrangement,  it  is  seen  that 


(2.11) 


which  is  a  result  that  is  identical  to  eauation  t2.3).  Of 
course  a  similar  development  will  provide  a  sxibstantiation 
of  equations  (2.4),  (2.5)  and  (2.6). 

2 .  Node  to  Datxim  Analysis 

Using  the  approximations  for  the  second  derivative, 
a  node  analysis  can  be  performed  to  predict  the  temperature 
distribution  in  a  mathematical  model  of  a  physical  configura¬ 
tion.  With  a  node  analysis,  the  configuration  of  interest 
is  divided  into  small  but  finite  subvoliames  with  each  sub¬ 
volume  being  considered  as  isothermal.  It  may  be  then 
assumed  that  the  geometric  center  of  each  siibvolume  is 
representative  of  that  subvolume  and  this  geometric  center 
is  called  the  node-point  or  node. 

A  node  analysis  consists  of  an  application  of  the 
First  Law  of  Thermodynamics  to  generate  energy  balances  for 
each  node.  This  is  done  for  each  node  by  equating 

HEAT  IN  +  HEAT  GENEEATED  =  HEAT  OUT  +  HEAT  ACCUMULATED  (2.12) 


For  the  two  dimensional  case  displayed  in  Figure  2.2,  the 
development  of  the  node  equations  for  node  point  (j,k)  is 
shown  in  Figure  2.3.  Assuming  heat  flow  by  conduction  and 
a  nodal  dissipation,  q,  a  steady  state  or  First  Law  energy 
balance  requires  that  the  Siam  of  all  heat  flows  (Watts  or 
BTU/hr)  leaving  the  node  equals  zero.  Thus  with  a  heat  input 
as  a  negative  heat  flow  outward, 

+  ^2  ^3  '*■  ^4  ~  ^^i  ~  °  (2.13) 

In  terms  of  the  temperatures  in  Figure  2.4  and  in  terms  of  the 
Fourier  Law 


KA(AT) 

^  '  T - 


(2.14) 


q 


2 


KAx 

Ay 


5  (T  -T,) 

O  2 


^3 


<54 


KAx 

Ay 


and  when  these  are  inserted  into  equation  (2.13),  the  result 


(2.15) 


KAy 

AX 


6 (T  -T, ) 
O  1 


KAx 

Ay 


6 (T  -T.) 
O  2 


6(T^-T3) 


+  6  (T  -T.)  =  q. 

Ay  o  4' 


If  two  conductances  are  defined  by  Kl  =  KAy6/Ax  and 
K2  =  KAx6/Ay  then 

KI(Tq-T^)  +  K2(Tq-T2)  +  K1(Tq-T3)  +  K2(T^-T4)  = 

(2.16) 

and  if  there  is  no  heat  generation  (q^=0) ,  the  equation 
becomes 


Tq(2K1  +  2K2)  -  Kl(T^)  -  K2(T2)  -  Kl (T3)  -  K2(T4)  =  0 

(2.17) 


For  the  case  where  there  is  heat  flow  by  both 
conduction  and  convection  as  in  Figure  2.5,  with  a  heat 
input  as  a  negative  flow  outward, 

q^  +  q2  +  q3  +  -  qi  =  0  (2.18) 

For  the  conduction  terms,  q^,  q2  and  q3,  equation  (2.14)  is 
still  applicable.  However,  for  the  convection  at  the  boundary 
as  indicated  by  q^  (Figure  2.5) 


hAy6 (AT) 


(2.19) 


In  terms  of  the  temperatures  in  Figure  2.5 


6  (T  -T,  ) 
Ly  o  1 


12  = 

,3  =  5  ,t^-T3) 

=  hAy  6 (T^-T^) 

and  when  these  are  inserted  into  equation  (2.18).,  the  result 
is 


+  hAy6 (T^-Tg)  =  q^  (2.20) 

Letting  K1  =  K2  =  and  K3  =  hAy6 ,  then 

KI(Tq-T^)  +  K2(T^-T2)  +  KKT^-T^)  +  K^(T^-Tg)  =  q^ 

(2.21) 

which  can  be  simplified  to 

(2K1  +  K2  +  K3)T  -  KlT  -  K2T-  -  KlT,  =  q.  (2.22) 


and  once  again  if  there  is  no  heat  generation  (q^  =0),  the 
equation  becomes 


(2K1  +  K2  +  K3)T^  -  KlT^^  -  -  KlT^  "  KST^  =  0  (2.23) 

For  the  case  of  a  node  on  an  adiabatic  boundary  with  no 
heat  generation  as  in  Figure  2.6 


:-v\ 

■\vtvv 


+  ^2  ^3  "  ° 


(2.24) 


In  terms  of  the  temperature  in  Figure  2.6  and  in  terms  of  the 


Fourier  Law 


6(T  -T-) 

Ax  ^  o  ^2' 

gI^x/2-)_  j 

Ay  o  3' 


and  when  these  are  inserted  into  equation  (2.24),  the  result 


"‘‘y  ^  «<To-T2>  *  = 


(2.25) 


.  ..  .  K(Ax/2)  r  -  2.  KAy6 

Letting  K1  =  — ^ — Ay^  °  “  ~Ax 


(2K1  +  K2)  Tq  -  KlTj^  -  K2T2  -  KIT^  =  0  (2.26) 


It  is  easy  to  see  that  when  one  proceeds  in  this 
manner  for  all  nodes,  an  nxn  set  of  simultaneous  algebraic 
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equations  will  result.  These  may  be  represented  in  matrix 
form  as 

[K] [T]  =  IQ]  (2.27) 

where  [K]  is  an  nxn  matrix  of  coefficients  and  IT]  and  IQ] 
are  nxl  column  vectors  representing  the  node  temperatures 
and  heat  inputs  respectively. 

The  solution  of  equation  (2.27)  is  handled  by  the 
computer  program  TVSSI  [Ref.  12; Appendix  A] . 

B.  MODEL  BUILDER  TVIN 

In  order  to  provide  input  to  program  TVSSI,  it  was 
necessary  to  create  a  model  of  the  configuration  to  be 
analyzed.  The  overall  model  of  the  condensing  region 
incorporated  three  major  features.  Because  the  system  to 
be  modelled  was  an  externally  finned  tube  with  the  external 
space  flooded  by  condensate,  the  azimuthal  portion  of  the 
model  was  chosen  to  be  l/27r  radians.  This  is  shown  in 
Figure  2.7  and  because  the  areas  of  each  node  contain  a  tt 
term,  the  use  of  l/27r  simplified  the  construction  of  the  node 
equations  by  peirmitting  cancellation  of  the  it  terms.  The 
second  feature  required  that  the  maximum  number  of  nodes  to 
be  used  for  the  model  was  300,  due  to  input  restrictions 
(the  dimension  statement)  of  the  main  thermal  analyzer.  The 
final  requirement  necessitated  that  the  program  be  written  so 


Tube  Wall 


Figure  2.7. 


Model  Azimuthal  Section 


as  to  allow  the  user  to  input  the  variables  listed  in  Table  I 
interactively . 

The  basic  model  is  shown  in  Figure  2.8.  The  upper 
horizontal  edge  represents  the  condensate/vapor  interface 
while  the  lower  horizontal  edge  represents  the  inner  wall  of 
the  tube.  Heat-transfer  coeffients  were  specified  for  the 
inner-tube  wall  and  the  vapor-condensate  interface.  The 
model  assumes  that  the  heat  is  transferred  across  the  retained 
condensate  by  conduction  only.  Due  to  symmetry,  the  dotted- 
dashed  lines  at  the  two  vertical  edges  of  the  figure  represent 
an  isolation  of  a  fin-condensate-tube  entity  into  a  repeating 
section.  These  lines,  of  course,  represent  adiabatic 
boundaries.  Thus,  the  model  assumes  that  the  only  heat 
transferred  into  or  out  of  the  model  occurs  at  either  the 
condensate/vapor  interface  or  the  inner  tube  wall. 

The  first  model  constructed  consisted  of  299  nodes  with 
13  columns  and  23  rows.  The  outermost  rows  and  columns 
consisted  of  half  node  subvolumes.  This  facilitated  the 
placement  of  the  node  centers  at  both  the  inner  tube  wall 
and  upper  condensate/vapor  interface  in  addition  to  the 
boundary  of  the  outer  tube  wall  and  fin  perimeter.  Nine  rows 
were  placed  in  the  condensate  region  above  the  fin  and  in  the 
region  including  the  fin  and  adjacent  condensate.  The 
remaining  four  rows  were  placed  in  the  tube  wall  region. 

In  addition  the  fin  region  contained  4  columns  of  nodes 
while  the  remaining  eight  columns  were  placed  in  the  inter-fin 
space . 


TABLE  I 


The  location  of  these  columns  and  rows  determined  the 


height  and  width  increments  of  the  various  subvolumes, 
depending  on  their  location  in  the  model.  For  example,  if 
the  tube  wall  thickness  is  four  millimeters,  and  the  tube 
wall  area  contains  four  rows,  the  height  increment  for  a 
node  in  this  region  would  be  4  mm/4  rows  =  1  mm.  The  width 
increments  for  a  node  would  be  computed  in  identical  fashion 
by  dividing  the  size  of  the  particular  region  by  the  number 
of  columns  occupying  that  region.  Figure  2.9  identifies  all 
of  the  different  height  and  width  increments  utilized  by  the 
model . 

C.  EVALUATION  OF  RESULTS 

1 .  Verification  of  Accuracy 

Verification  of  the  model  consisted  of  using  the 
temperatures  generated  by  TVSSl  (Using  TVIN  as  input)  and 
the  conductance  values  employed  in  TVIN  to  compute  the  heat 
flow  across  three  boundaries  in  the  model : 

(1)  condensate  upper  surface 

(2)  fin  perimeter  and  tube  outer  wall 

(3)  tube  inner  wall 

The  foregoing  heat  flows  were  computed  in  each  case 
by  using  the  node  temperatures  generated  by  TVSSl  with  the 
associated  internode  conductance  values  used  in  the  input 
generated  by  TVIN.  These  heat  flows,  of  course,  must  be 


Height  Increment  3 


Condensote 
Upper  Level 


Width  Increment  I  Width  Increment  2 


B 

Figure  2.9.  Model  Increments.  (A)  Height  and  (B)  V7idth 
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equal  and  it  was  observed  at  the  outset  that  13x23  model 
yielded  certain  unacceptcible  discrepancies  in  the  three  values 

Because  the  13x23,  299  node  model  was  found  to  be 
inaccurate,  it  was  determined  that,  in  order  to  achieve  the 
requisite  accuracy,  slight  adjustment  of  the  number  of  rows 
and  columns  was  required  from  run-to-run. 

TVIN  was  reprogramed  to  allow  flexibility  in  the 
number  of  rows  and  columns  utilized  in  the  model  as  well  as 
placement  of  the  nodes  in  the  fin,  tube  or  liquid  regions. 

The  first  change  instituted  was  to  fix  the  number  of  rows 
contained  in  the  condensate  layer  above  the  fin  tip  at  two 
rows.  This  action  was  acceptable  because  the  original  height 
of  this  layer  was  only  0.08  mm  compared  with  a  1.0  mm  fin 
height  and  3.175  mm  tube  wall  thickness.  Fixing  the  number 
of  rows  at  two  permitted  the  use  of  additional  rows  at  the 
fin  and  tube  levels.  Program  TVIN  was  changed  to  inter¬ 
actively  accept  a  specified  number  of  columns  (see  Appendix 
A) .  The  number  of  rows  was  obtained  by  making  the  product 
of  rows  and  columns  as  close  as  possible  to  300. 

A  description  of  the  menu-driven  procedure  for  the 
selection  and  allocation  of  rows  and  columns  for  a  296-node 
model  example  now  follows: 

If  the  response  to  the  query  'How  many  columns  are 
desired?'  is  8,  the  program  will  provide  300/8  =  37.5  or  37 
rows.  This  yields  an  8x37  =  296  node  model. 

Suppose  the  response  to  the  query  'How  many  rows  are 
desired  for  the  tube  wall?'  is  5.  Then  the  program  will 


allocate  37-5  =  32  rows  for  the  fin  and  condensate  layer 
above  the  fin. 

However,  the  number  of  rows  in  the  condensate  layer 
above  the  fin  tip  is  fixed  at  2.  Thus  the  number  of  rows  in 
the  fin  becomes  32-2  =  30.  The  result  is  an  8x37,  296  node  model 
A  display  of  the  foregoing  model  set-up  is  shown  in 
Figure  2.10. 

2 .  Model  Verifier  Q1 

As  mentioned  in  the  previous  section,  once  an  input 
file  had  been  created  by  TVIN,  and  the  node  temperatures 
were  calculated  by  the  thermal  analyzer,  TVSSI,  the  solution 
was  verified  for  correctness  and  accuracy.  This  check 
consisted  merely  of  a  comparison  of  the  heat  flow  (watts) 
across  three  surfaces  in  the  model.  An  error  of  3.5%  was 
established  as  the  meiximum  allowable  for  a  solution  to  be 
considered  acceptable.  The  model  verifier  Program  Ql  (see 
Appendix  B) ,  is  an  offshoot  of  the  model  builder  TVIN  as  it 
interactively  accepts  the  same  input  data  required  to  build 
the  model.  Once  all  input  parameters  are  set,  the  program 
generates  and  stores  each  inter-connecting  node  conductance 
value . 

In  the  next  step,  all  node  temperatures  computed  by 
the  thermal  analyzer  are  read  and  stored,  by  location  in  the 
model.  The  node  numbers  of  the  condensate  outer  layer,  fin 
perimeter  and  outer  tube  wall  and  the  inner  tube  surface  are 
determined  for  each  model  via  an  algorithm  which  requires  as 


input,  the  number  of  columns  and  rows  being  used  in  the  model 
along  with  their  placement.  In  addition  to  determining  the 
nodes  on  these  surfaces,  the  algorithm  also  computes  the 
node  nximber  of  each  of  the  adjoining  nodes  in  the  direction 
of  the  heat  flow.  Computation  of  the  heat  flow  by  program  Q1 
is  then  simply  a  matter  of  taking  each  surface  node  conductance 
and  multiplying  it  by  the  temperature  difference  between  the 
surface  node  and  its  adjoining  node.  The  summation  of  all 
the  surface  node  values  results  in  the  total  heat  flow  across 
the  particular  surface. 

D.  THERMAL  ANALYZER  TVSSI 

The  computer  program  TVSSI  used  in  obtaining  the  temperature 
distributions  was  adapted  from  program  TVSS2  listed  by  Kern 
and  Kraus  [Ref.  13: Appendix  2a].  This  adaptation  consisted 
of  modifying  TVSS2  to  compute  solutions  in  the  SI  system  and 
accept  an  input  file  specifically  designed  by  model  builder 
TVIN  for  each  externally  finned  tube  case.  This  program  also 
had  to  be  adapted  to  the  IBM  3033  AP  system. 

The  program  utilizes  the  Cholesky  decomposition  scheme, 
and  because  of  the  linearization  of  the  radiation  terms  (not 
used  in  this  study),  the  program  is  iterative.  A  more  indepth 
description  of  the  Cholesky  Scheme  pertaining  to  the  program 
is  discussed  by  O'Hare  [Ref.  12:pp.  48-52]. 

For  all  the  cases  which  were  studied,  the  following 
parameters  were  set  for  the  analyzer: 


1.  The  initial  temperatures  at  which  all  nodes  were  set 
and  at  which  the  computer  began  the  iterative  process 
was  100°C. 

2.  An  eventual  accuracy  of  .05  between  the  final  and  next 
to  last  iterations  was  used. 

3.  A  radiation  coefficient  convergence  factor  Of  0.66667 
between  iterations  was  used.  This  is  required  by  the 
program  even  if  radiation  is  not  considered  as  a  mode 
of  heat  transfer  in  the  model. 

4.  The  maximum  number  of  iterations  that  the  computer  was 
allowed  to  perform  was  set  at  12  to  avoid  excessive 
computer  running  time. 

5.  A  damping  factor  of  0.8  was  set  as  an  initial  dcimping 
factor  based  on  the  number  of  non-linear  terms  in  all 
of  the  node  equation. 

The  ability  of  TVSSI  to  correctly  solve  a  numerical  model 
was  verified  by  providing  as  input  a  data  set  listed  by  Kern 
and  Kraus  [Ref.  13;pp.  326-336]  for  an  extended  surface 
problem  where  the  solutions  were  derived  analytically.  The 
computer  generated  solution  agreed  completely  with  the 
analytic  solution  thereby  assuring  confidence  of  analyzer 
generated  temperatures. 


III.  RESULTS  AND  DISCUSSION 


A.  INTRODUCTION 

During  this  thesis  effort,  15  different  cases  were  studied 
using  the  programs  TVIN  and  TVSSI.  Table  II  lists  the  param¬ 
eters  that  remained  constant  in  all  cases.  A  sxammary  of  all 
cases  examined  and  their  results  are  listed  in  Table  III. 

B.  SMOOTH-TUBE  EVALUATION 

As  a  method  for  comparison  of  the  variation  in  heat- 
transfer  performance,  the  heat  flux  for  a  smooth  tube  with 
the  same  diameter  as  the  finned-tube  root  diameter  was  estab¬ 
lished  as  a  baseline  value.  The  ratio  of  the  one-  and  two- 
dimensional  heat  fluxes  to  the  heat  flux  of  the  finned  tube 
with  fin  width  approaching  zero  were  designated  and  (^2 
respectively.  The  equation  used  to  determine  the  finned-tube 
of  zero  fin  width  heat  flow  is  given  by: 

(Tq  -T^) 

0  =  D^(ln(DyD^))  ,  e  ,  1  ».l) 

(h^)  (D^)  2k„ 

C.  ONE  DIMENSIONAL  HEAT  FLOW 
The  one-dimensional  model  utilized  in  this  study  assumed, 

as  did  Owen  et  al.  [Ref.  8],  parallel  paths  for  heat  flow  across 
the  tube  wall,  fin  and  retained  condensate.  Furthermore,  a 
simplifying  assxamption  was  made  that  heat  flow  is  purely  radial 
and  that  the  heat  is  transferred  across  the  retained  condensate 


TABLE  II 


I 

I 


CONSTANT  PARAMETERS  FOR  MODEL  BUILDING  PROGRAM  TVIN 


h^  =  30000  W/m-K 

h  =  100000  W/m-K 

o 

=  6 . 35  mm 

r  =  9.525  mm 

o 

e  =  1 . 0  mm 


K^  =  0.6  W/m-K 

T  =  48.0  "C 

o 

T^  =  22.0  ”€ 
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SUMMARY  OF  CONP IGURATIUNS  EXAi>lINLU  WITu  HEAT  TRANSFER  PERFORl-J\NCE 
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by  conduction  only.  This  was  accomplished  by  modifying  pro¬ 
gram  TVIN  in  such  a  way  so  that  the  data  file  it  created  con¬ 
tained  only  those  internode  conductance  values  representing 
heat  flow  in  the  radial  direction  while  deleting  those  conduc¬ 
tance  terms  pertaining  to  axial  heat  flow.  Programs  TVSSI  and 
Q1  were  utilized  to  compute  the  node  temperatures  and  heat 
flow  for  each  case.  Once  again,  a  maximum  error  of  3.5%  was 
used  when  heat  flows  generated  by  program  Ql  were  verified. 

D.  TWO  DIMENSIONAL  HEAT  FLOW 

The  heat  flow  for  the  two-dimensional  case  was  evaluated  by 
the  methods  discussed  in  Chapters  I  and  II. 

E.  EFFECTS  OF  FIN  SPACING  AND  FIN  WIDTH  ON  PERFORMANCE 

Heat  fluxes  were  computed  for  five  different  cases  of  fin 
width  and  spacing.  Three  cases  had  a  fin  width  of  1  mm  with 
spacings  of  0.5,  1.0  and  1.5  mm  respectively.  The  other  two 
cases  held  the  spacing  constant  at  0.5  mm,  while  the  fin  width 
was  adjusted  to  0.5  and  0.25  mm. 

1 .  Effects  of  Fin  Spacing 

This  section  presents  results  showing  the  variations 
of  heat  flux  having  fin  spacing  as  a  parameter.  Figures  3.1, 

3.2  and  3.3  present  results  for  tubes  having  a  fin  width  of  1.0 
mm  with  fin  spacings  of  0.5,  1.0  and  1.5  mm.  Each  figure  repre¬ 
sents  the  same  spacing  variation  with  the  tube-metal  thermal 
conductivity  being  varied  from  385.0  W/m.K  for  copper,  through 
45.0  W/m  K  for  copper-nickel,  down  to  16.0  W/m-K  for  the  case  of 
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Predicted  Heat-Transfer  Performance  of  Finndd  Copper 


Figure  3.2.  Predicted  Heat-Transfer  Performance  of  Finned  Copper 
Nickel  Tube  With  e  =  1.0  mm  and  t  =  1.0  mm 


Figure  3.3.  Predicted 
Steel  Tub 


stainless  steel.  Similarly,  Figures  3.4,  3.5  and  3.6  present 
the  ratio  of  the  two-dimensional  heat  flux  to  the  one-dimensional 
heat  flux  (Q21) .  For  each  case  of  tube-metal  thermal  conduc¬ 
tivity,  the  fin  spacing  of  0.5  mm  shows  the  best  heat-transfer 
performance,  while  the  worst  performance  occurs  at  a  spacing  of 
1.5  mm.  Figures  3.7,  3.8  and  3.9  present  isotherm  plots  for 
all  three  cases  of  thermal  conductivity  for  the  finned  tube 
have  a  0.5  mm  fin  spacing  and  1.0  mm  fin  width.  Moreover,  the 
ratio  Q21  has  its  smallest  value  at  a  find  spacing  of  0.5  mm 
and  increases  with  increasing  fin  spacing.  Enhancement  ratios 
for  the  various  cases  had  their  maximum  value  at  the  0.5  mm 
fin  spacing  and  this  ratio  showed  a  marked  decreaase  as  the  fin 
spacing  was  increased. 

The  results  indicate  several  important  facts.  The  first 
and  most  important  is  the  significant  role  the  two-dimensional 
heat  flow  plays  in  this  type  of  extended  surface  analysis.  The 
increasing  Q21  ratio  with  increasing  fin  spacing  shows  that  for 
the  cases  of  copper  and  copper-nickel  tube-metal,  the  axial 
heat  flow  into  the  fin  upper  side  surface  tends  to  raise  the 

temperature  of  the  fin  tip  resulting  in  a  choking  effect  of  the 

heat  flow  into  the  fin  tip  surface.  This  is  of  enormous  conse¬ 
quence  because  in  each  and  every  case,  the  majority  of  the 

he  -  absorbed  by  the  finned  tube  is  accomplished  at  the  fin 
tip.  Because  the  one-dimensional  case  allows  no  axial  heat 
flow  into  the  fin  side  surface,  a  greater  temperature  gradient 
occurs  at  the  fin  tip  resulting  in  a  greater  heat  flow.  As  the 
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Figure  3.6.  Comparison  of  One-  and  Two-Dimensional  Heat-Transfer  Perfor¬ 
mance  For  Stainless  Steel  Tube  With  e  =  1.0  mm  and  t  =  1.0  m 


Figure  3.8 


Isotherm  Plot  For  Finned  Copper-Nickel  Tube 
(Fin  and  Condensate  Regions) 


Fiaure  3.9. 


Isotherm  Plot  For  Finned  Stainless  Steel  Tube 
(Fin  and  Condensate  Regions) 


fin  spacing  is  increased,  the  effect  of  the  fin  side  surface 
is  less  pronounced  since  the  poor  conductivity  of  the  condensate 
inhibits  axial  heat  flow  in  the  region  near  the  fin  space 
centerline.  This  results  in  a  more  closer  agreement  between 
the  one-dimensional  model  and  the  two-dimensional  version. 

This  conclusion  can  be  further  justified  by  comparing 
the  Q21  ratios  for  each  tube-metal  conductivity.  The  two- 
dimensional  effect  at  the  fin  tip  is  decreased  as  the  thermal 
conductivity  of  the  fin  material  decreases.  This  does  not 
mean,  however,  that  there  is  no  two-dimensional  effect  for  the 
cases  with  low  tube-metal  conductivity.  As  in  the  case  of 
copper,  the  two-dimensional  effect  tends  to  inhibit  the  heat 
flow  into  the  fin.  For  the  cases  of  low  metal  thermal  conduc¬ 
tivity,  there  is  a  less  pronounced  two-dimensional  effect  at 
the  fin  tip,  but  there  is  a  significant  two-dimensional  heat 
flow  into  the  lower  portion  of  the  fin  side  surface  near  the 
tube  wall.  For  this  reason,  the  one-dimensional  model  tends 
to  underpredict  the  two-dimensional  model. 

As  seen  in  all  the  cases  for  the  copper-nickel  and 
stainless-steel  tube-metal  (except  the  copper-nickel  case 
with  0.5  mm  fin  spacing),  all  Q21  ratios  are  greater  than 
unity.  For  the  copper-nickel  case  of  0.5  mm  spacing,  the  reduced 
spacing  even  with  a  lower  metal  conductivity  of  45  W/m.K,  tends 
to  produce  a  two-dimensional  effect  at  the  upper  fin  tip. 

2 .  Effects  of  Fin  Width 

The  results  of  the  variation  of  heat  flux  with  fin 
width  as  a  parameter  indicate  that  fin  width  plays  an  important 


role  in  the  overall  heot- transfer  mechanism.  The  two  cases 
for  each  tube-metal  thermal  conductivity  investigated  were  fin 
widths  of  0.5  and  0.25  mm  respectively  with  fin  spacing  held 
constant  at  0.5  mm.  For  each  case  of  tube-metal  thermal  con¬ 
ductivity,  the  fin  width  of  0.5  mm  displayed  one  of  the  highest 
heat  fluxes,  being  second  only  to  the  case  of  the  1.0  mm  fin 
width  with  0.5  mm  fin  spacing.  The  only  significant  two- 
dimensional  fin  tip  effect  can  only  be  observed  in  the  case 
of  the  copper  tube-metal  with  a  fin  width  of  0.5  ram. 

The  same  reasoning  as  in  the  previous  section  can  be 
applied  in  this  case.  Because  the  fin  tip  width  has  been 
reduced,  the  total  amount  of  heat  transferred  into  the  tip 
region  is  also  reduced  and  the  two-dimensional  fin  tip  choking 
effect  has  less  influence.  This  concurs  with  the  computed  Q21 
ratio  of  1.04.  By  reducing  the  fin  width  further  in  each  case 
to  0.25  mm,  each  of  the  resultant  heat  fluxes  is  less  than  the 
cases  with  0.5  mm  fin  width.  For  each  case  of  the  0.25  mm  fin 
width  the  Q21  ratio  approached  1.3,  signifying  the  comparatively 
large  two-dimensional  heat  flow  into  the  lower  fin  side  surface. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  The  two-dimensional  niomerical  model  of  a  horizontal 
finned  txibe  fully  flooded  with  steam  condensate  developed  is 
a  useful  tool  in  the  study  of  the  relationship  between  one- 
and  two-dimensional  heat- transfer  effects. 

2.  Two-dimensional  heat  transfer  effects  play  a  major 
role  in  the  overall  heat  transfer  mechanism  in  the  case  of  an 
externally  finned  tube  fully  flooded  with  steam  condensate. 

It  has  been  determined  that  for  cases  of  high  tube-metal  con¬ 
ductivity,  such  as  the  case  of  copper,  the  one-dimensional 
model  can  overpredict  the  results  of  the  two-dimensional 
model  by  as  much  as  13%,  while  variation  in  the  fin  width  or 
spacing  can  result  in  the  one-dimensional  model  underpredicting 
the  two-dimensional  result  by  as  much  as  twenty-two  percent. 

3.  Fin  width  and  fin  spacing  control  the  magnitude  of  the 
two-dimensional  effect  present  and  the  resulting  heat- transfer 
performance.  An  increasing  fin  width  to  fin  spacing  ratio 
will  correspond  to  an  increase  in  the  heat- transfer  performance. 

4.  The  one-dimensional  model  in  cases  of  high  tube-metal 
thermal  conductivity  exhibiting  a  fin  width  to  spacing  ratio 

of  greater  than  unity  will  overpredict  the  heat  flux  predictions 
of  the  two-dimensional  model. 

5.  Decreasing  the  tube-metal  thermal  conductivity  results 
in  a  degraded  heat-transfer  performance. 


6.  The  two-dimensional  fin-tip  choking  effect  is  minimal 
in  low  thermal  conductivity  tube-metals  while  the  two- 
dimensional  heat  transfer  into  the  fin  lower  side  surface  is 
far  more  pronounced. 

B.  RECOMMENDATIONS 

1.  Use  the  two-dimensional  model  computed  enhancement  ratios 
for  the  various  cases  to  estimate  a  steam-side  heat-transfer 
coefficient  and  compare  with  those  methods  currently  being  used. 

2.  Using  the  model  builder  TVIN  and  thermal  analyzer 
TVSSI,  conduct  a  more  in-depth  study  of  the  effects  of  variation 
of  fin  width,  height  and  spacing  on  the  overall  heat-transfer 
performance. 

3.  Revise  the  two-dimensional  model  to  include  convection 
effects  in  the  retained  condensate  region. 
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TVIN  COMPUTER  GENERATED  INPUT  ANALTZER  PROGRAM 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM  TVIN 

PROGRAM  TVIN  IS  DESIGNED  TO  ACCEPT  VARIABLE 
DIMENSIONS  AND  INPUT  PARAMTERS  TO  CREATE 
AN  OUTPUT  DATA  FILE  OF  INTER  NODE  CONNECTIONS 
AND  CCNOUCTANCE  TERMS.  BY  ADJUSTING  THE 
PRINT  STATEMENT,  THE  OUTPUT  WILL  EITHER 
REPRESENT  A  ONE  OR  TWO-DIMENSIONAL  MODEL. 


CCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


DIMENSION  CK(300,4)tT(300) 

INTEGER  P(300.4)  ,C0,N,X,ZW,ZW1,ZW2,ZW3,V1,V2 ,W,Uf 
DOUBLE  PRECISION  CTHM ,CTHL, HI ,HS , RO , RI , A ,B, C , 0, 
♦VAL4,VAL14,VAL6,VAL16.VAL8,VAL18,VAL5,VAL7, VALID, 
*A2,A3iA4,A5,WINCl , WINC2 , H INCl ,HI NC2, HI NC3 , AL INE , 
4<TIN,  TS,  T,X1  ,X2,  VAL2,VAL12,VAL20,  AI  ,RADH,RADL  ,EL 
CHARACTER  *79  TITLE 
CHARACTER  *79  FNAME 
C 

DATA  IOI,IN,  IPR,IWR/6,5,4,8/ 

C 

OP£N( IPR,FILE='PRN',STATUS='NEW'  , FQRM  =  • FORMATTED* , 
*IOSTAT=ICK} 

IFdCK  .NE.  0)  WRITEl  lOT  ,20  ) 

^20  FORMAT! 'TROUBLE  OPENING  PRINTER  OUTPUT  FILE*) 

C 

WRITEilOT, 17) 

17  FORMAT!/'  ENTER  THE  TITLE  OF  THIS  STUDY,  79  *, 

♦  *  COLUMNS  MAX.  THIS  TITLE  WILL  APPEAR*,/,*  AI  * , 

♦  •  THE  TOP  OF  EVERY  PRINTED  PAGE  OF  OUTPUT:*) 
REA0!IN,18)  TITLE 

18  F0BMAT!A79) 

C 

WRITEIIOT, 1 ) 

1  FORMAT!/*  ENTER  INNER  RADIUS  OF  TUBE:') 

REA0!IN,2)  RI 

2  F0RMAT!BN,F10.5) 

C 

WR1TE1I0T,3) 

3  FORMAT!/'  ENTER  OUTER  RADIUS  OF  TUBE:*) 

REA0!IN,2)  RO 

C 

WRITE!I0T,4) 

4  FORMAT!/*  ENTER  DISTANCE  FROM  OUTER  EDGE  OF  *, 

♦  'FIN  TO  LIQUID  SURFACE:*) 

REA0!IN,2)  A 

C 

rtRITE!I0T,5) 

5  FORMAT!/*  ENTER  HEIGHT  OF  FIN*) 

READ!IN,2)  B 

C 

WRITE!ICT,6) 

6  FORMAT!/*  ENTER  DISTANCE  BETWEEN  FINS') 

READ1IN,2)  C 
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WRITEdOT,?) 

FORMATC/'  ENTER  WIDTH  CF  FIN*) 

READ(IN,2)  0 

WRITEdOT, 8) 

FORMAT!/'  ENTER  HEAT  TRANSFER  COEFFICIENT  OF 
*•  TUBE  INNER  SURFACE:  •  ) 

READ(IN,2)  HI 

WRITEdOT, 9) 

FORMAT!/'  ENTER  HEAT  TRANSFER  COEFFICIENT  OF  ', 

♦  'OUTER  LiaUIQ  SURFACE:*) 

REAO!IN,2)  HS 

WRITE!I0T,10) 

FORMAT!/'  ENTER  THERMAL  CONDUCTIVITY  OF  FINS  AND  ', 

♦  '  TUBE:') 

READ(IN,2)  CTHM 

WRITEdOT, 11) 

FORMAT!/'  ENTER  THERMAL  CONDUCTIVITY  OF  FLUID  '  , 

♦  '  SURROUNDING  TUBE  AND  FINS:') 

READ!IN,2)  CTHL 

WR1TE!I0T,12) 

FORMAT!/'  ENTER  TEMPERATURE  OF  FLUID  INSIDE  TUBE:') 
READIIN,2)  TIN 

WRITEdOT,  13) 

FORMAT!/'  ENTER  TEMPERATURE  EXTERNAL  TO  *, 

*'  SURROUNDING  FLUID:') 

READ! IN, 2)  TS 
WRITE!I0T,14) 

FORMAT!/'  ENTER  THE  NUMBER  OF  COLUMNS  DESIRED  ', 

FOR  YOUR  model:  ') 

READdN,15)  N 
F0RKAT!BN,I2) 

WRITEIIOT,  16) 

FORMAT!/*  ENTER  THE  NUMBER  OF  ROWS  DESIRED  FOR  ', 

*•  THE  TUBE:') 

READ! IN,15)U 
WRITEdOT, 99) 

FORMAT!/'  ENTER  THE  NUMBER  OF  COLUMNS  DESIRED  FOR  ', 
*'  THE  fin:  ') 

REAO!IN,15)ZW 

C0=INT!300.0/FL0AT!N)  ) 

INITIALIZE  THE  P  AND  CK  ARRAYS 

DO  25  I=l,CO*N 
DO  25  J  =  l,4 
P!  I,  J)  =  0 
CK  !I  ,U)=0.0 
CONTINUE 

VALUES  COMMON  TO  ALL  RUNS 

L21=CO*N 

L22=2 

L23=0 


v'f' 

V'l-'v 


■:» 
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L24»0 

L25=0 

L26=0 

L27=0 

L28=0 

L31=300 

L32=50 

L33=6 

L34=2 

L3  5=4 

L36=6 

FL41=.01 
FL42=. 66667 
IL43  =  12 
FL44  =  .8 
FL45=100. 

LI  NT  1=0 
LINT2=0 
LI  NT  3=0 

COMPUTE  INTERNODE  CONNECTIONS 

DO  55  I  =  1,C0*N 
X  =  N*(C0-1)  +  1 

IF(AMOO(FLOAT( I-l I .FLOATtN) I.EQ.O)  GO  TO  30 
IF|AMODIFLOAT( I),FLOAT(N) J.EQ.O)  G3  TO  40 
IFil.LE.N)  63  TO  42 
IF(l.GE.X)  60  TO  48 


PIIfl)  =  U-N)4lO+l 
P(  I,2J  =  (I-IJ*10+1 
Pdf  3)=(I  +  1)*10+1 
P(It4)=(I+N)*10+l 
GO  TO  55 


IF(I.LE.N)  63  TO  45 
IFd.GE.XJ  60  TO  49 
P(  I,2)=(  I  +  U*10+l 
P(  If  i)=(  I-N)*10+l 
P( I,3i=(I+N)*lO+l 
GO  TO  55 


IF(I.GE.X)  GO  TO  51 
IF(I.LE.N)  63  TO  47 
P4I,2)=(I-1I*10+1 
GO  TO  35 


P(Ifl)={I-l)*10+l 
P  ( If  2)=d  +  ll*10+l 
P(I,3)=(I+N)*10+1 
P( If4)=3014 
GO  TO  55 


P(  If  li={ I*20)  +  l 
P<  If 2)=l I  +  N)*10+l 
P(  If 3)=3014 
GO  TO  55 


GO  TO  150 


IF(l.EQ.CO)  GO  TO  125 

RAOH=RO*B+(fLOAT( I )-( FLO ATC U) +0. 5 ) -FLOAT ( ZW3 ) INC3 
RAOL=RO-frB+(FLOAT(  I  J-(FL3AT{  U)  +  l.  5)-FL0AT  (  ZW3  )  )*HINC3 
GO  TO  150 

RAOH=RO+B+A 

RA0L=R0+fa+A-HlNC3/2.0 


AL INE=RADH**2-RA0L**2 

COMPUTE  INTER-NOOE  CONDUCTANCE  VALUES 


IF(I.GT.l)  GO  TO  181 

HTC=HI 

RD=RI 

GO  TO  189 

iFCI.NE.COI  GO  TO  1< 

HTC=HS 

RD=RO+A+E 

Al=HTC*WINCl*R0/2 .0 

A2=Al*2-0 

A3=A2*eL/WINCl 

A4=HTC*WINC2*RD 

A5=A4/2.0 


IF(I .GE. (U+2J )  GO  TO  203 

VAL4=RADL*WINC1*CTHM/H1NC1 
GO  TO  300 

IF<I.G£.(  (U+2)+Zw3n  GO  TO 
VAL4=RA0L*W1NC1*CTHM/HI'JC2 
GO  TO  300 

VAL4=CTHL*RADL*WI NC1/HINC3 
VAL2=VAL4/2.0 


IF(I  .G£.(U+1J )  GO  TO  310 
VAL14=CTHN*RA0H*W INCI/HINCI 
GO  TO  400 

IF(I .LE.(U+ZW3 )»  GO  TO  350 
VAL14=CTHL*RA0H*M INC1/HINC3 
GO  TO  40u 

VAL14=C THM*RA0H*WINC1/HINC2 
VAL12=VAL14/2.0 


1F(1  .GE.(U-«-2)  )  GO  TO  413 
VAL6=CTHM*RA0L*EL /HINCl 
GO  TO  50C 

IF(I  .GE.(<U+2J+ZW3i)  GO  TO 

val6=cTcthm*winci+cthl*wini 

GO  TO  500 

VAL6=CTHL*RADL*EL/HINC3 


THL*WINC2J*R40L)Z (2.0*HIMC2) 


70 


IF(I.GE.(U+in  60  TO  613 
VAL16=CTHM*EL*RA0H/HINC1 
GO  TO  650 

IF(I.LE.(U>ZM3]}  GO  TO  620 
VAL16=CTHL*EL*RADH/HINC3 
GO  TO  650 

VAL16=I  (CTHM*WINC  1+CTHL*W  INC2  »*RADH  )/  (2.0*HINC2  ) 


IF(l.GE.(U+2i J  60  TO  670 
VAL8=RA0L*WINC2*CTHM/HINCl 
GO  TO  700 

IF(I.GE.((U+2)>ZM3)I  GO  TO 
VAL8=CTHL*RA0L*WINC2/HINC2 
GO  TO  700 

VAL8  =CTHL*RA0L*WINC2/HINC3 
VAL10=VAL8/2. 0 


IF(I  .GE.IU+li)  GO  TO  730 
VAL18=CTHM*RADH’^VI  INC2/HINC1 
GO  TO  750 

IF(I.LE.(U+ZW3))  0  TO  740 
VAL18=CTHL*RADH*W1NC2/HINC3 
GO  TO  750 

VAL18=CTHL*RADH*W INC2/HINC2 
VAL20=VAH8/2.0 


IF  (I. N£.<(U+1) +711(3)1  GO  TO  760 
VAL5=<CTHL*(RA0H**2-(R0+B)**2)+CTHM*nR0+B) 
♦♦♦2-RA0L**2J)/(2.0*WINCl) 

GO  TO  800 

IF(I.GT.((U+1)+ZW3))  GO  TO  770 
VAL5=CTHM*ALINE/(  2.0*WIMCl) 

GO  TO  800 

VAL5=CTHL*ALINE/I  2.0+WINCl) 


1F(I .GE.(U+1) )  GO  TO  810 
VAL7=CTHM*ALINE/( 2.0*WINC2) 

GO  TO  845 

IFU.NE.IU+I)  )  GO  TO  820 

VAL7=<CTHL*(RADH**2-R0**2)+CTHM*(R0**2-RADL<‘*2)  ) 
♦  /  (2.0+WINC2 I 
GO  TO  845 

VAL7=CTHL*ALINE/(  2.0*WINC2) 

COUPLE  CONDUCTANCE  VALUES  TO  CONNECTION  TERMS 

IF(I.EQ.l)  GO  TO  905 
IF(I.EO.CO)  GO  TO  912 
0=N*(I-1) 

00  900  J=0+lt0+N 


Iff  I  .Eu.  (0  +  11 i  GO 
IFK'.LT.iQ+ZR+lJ) 
IF(o.EO.(W  +  Z«+I)) 
IFU.NE.  (0+N)  )  GO 


TO  875 
GO  TO  880 
GO  TO  885 
TO  890 


CK(J.i)  =  VAL13 
CK(J,2)=VAL7 
CK<J,3)=VAL20 
60  TO  900 


CK(J,IJ =VAL4 
CK(J,2)=VAL5 
CK(J,3I=VAL5 
CK(J#4J=VAL14 
GO  TO  900 


CK(  J,1J  =  VAL6 
CK(J»2)-VAi.5 
CK(Jt3i=VAL7 
CK(Jt4)=VAL16 
GO  TO  900 


CK(J,1)=VAL8 
CK(J«2)=VAL7 
CK(  J,3i  =  VAL7 
CK(J,4)=VALL8 
GO  TO  900 


CK(J,IJ=VAL2 

CK(J,2)=VAL5 

CK(J,3)=VAL12 

CONTINUE 
GO  TO  1000 


DO  910  J=l,N 
IF(J.GT.l)  GO 
CK(  J,1J=VAL5 
CK(J,2)=VAL12 
CK<J,3)=A1 
GO  TO  910 


IF(J.GE.ZW+1J 
CK(J,1)=VAL5 
CK(J,2)=VAL5 
CK( J,31=VAL14 
CK(Jf4i=A2 
GO  TO  910 


TO  907 


IF(J.NE.(ZW+in  GO  TO  908 

CK( J,1)=VAL5 

CK(J,2J=VAL7 

CK(Jt3}=VAL16 

CK( J,4I=A3 

GO  TO  910 


IF(J.EQ.N)  GO 
CK(J,1|=VAL7 
CK( J,2)=VAL7 
CK(Jt3i=VAL18 
CK(  J,4)=A4 


TO  909 


i.-i  J.«u^  «» 


917 

1000 

C 

c 

c 


GO  TO  910 


CK( J.1J=VAL7 
CK(Jt2)=VAL20 
CK(  J»3)=A5 

CONTINUE 


GO  TO  1000 


Q=(C0-1J*N 
00  917  J=Q+1,Q+N 
IF(J.GT.(Q+li  )  GO 
CKiJf  1)=‘VAL2 
CK( Jf2)=VAL5 
CK(J.3}=A1 
GO  TO  917 


1F<J  .GE.(ZW4-a  +  l)) 
CK( J,1]=VAL4 
CK(J,2»=VAL5 
CKU.3)=VAL5 
CK(J.4)=A2 
GO  TO  917 


IF(J .NE. (ZM+Q+1) ) 

CK<J«11::VAL6 

CK(J«2)-VAL5 

CK<J,3J=VAL7 

CK(J«4)=A3 

GO  TO  917 


1F(J.EQ.(Q+N}  ) 

CK(J,li=VALd 

CK(J,2)=VAL7 

CK(Jt31=VAL7 

CK(J,4J=A4 

GO  TO  917 


CK(J«li=VAL10 
CK(J,2i=VAL7 
CK( J,31=A5 
CONTINUE 
CONTINUE 


TO  914 


GO  TO 


END  OF  CONDUCTANCE  CALCULATIONS 
WRITE(I0T,922  ) 

FORMAT!/, ^  ENTER  NAME  OF  INPUT  FILE':) 

REA0(IN.923  )FNAME 
F0RMAT(A12) 

CREATE  DATA  FILE  FOR  TVSSI 

OPEN(lWR,FILE=FNAME,STATUS='NEW'  ,FOR  M=  •  F  ORMATTE  O'  , 
♦  IOSTAT=ICK) 


c 


URITE(I0T,925) 

FORMAK/*  WRITING  FILE  TVSSI*) 


WRITE(I0T«926)C0*NiNfCQ.ZW2 

FORMAT!/*  THIS  'tlSf*  NODE  MODEL  HAS  *,I2f'  RO/!S  • 


ROMS* ) 


ROWS 


**  AND  *,I2f*  COLUMNS.  THERE  ARE  •,I2»'  ROWS*: 
WR1TE(I0T«927)ZW3 

FORMAT!*  IN  THE  CONDESATE  LAYER  AND  *,I2,*  ROl 
**  IN  THE  FIN.  THE  FIN  CONTAINS') 

MRITE!I0T«928)  ZW 
FORMAT!  12.*  COLUMNS.  *) 

WRITE  STATEMENTS 

WRITE!IUR«940)  TITLE 

WRITE !IWR. 944) L21 , L22«L23tL 24tL2 5. L26 . L27» L28 
WRITEIIUR.950 J  LINT1.LINT2.LINT3 
MRITE!IWR.  946)L31,L32.L33»L34.L35»L36 
WRITE !I WR, 947 ) FL4 1 .FL42, 1L43.FL44. FL45 
WR1TE!IWR.943)TI.  TS 


10  05 


1011 


1015 

C 

1100 

C 

C  f 

c 

94  0 
941 
94  2 

943 
C 

944 

945 

946 
94  7 

948 

949 

950 
960 
980 
C 

C 


DO  1100  i= 
If  (P(I.4J 
IFICKd  .1) 
IF!CK!I,3) 
WR1TE!IWR« 
MRITEIIWR. 
GO  TO  1100 
WRITE!1WR« 
WRITE!IWR, 
GO  TO  1100 
IFICKII.l) 
IF!CK!I,3) 
WRITE!IWR» 
WRITE!IWR, 
GO  TO  1100 
MRITE!IWR« 
MRITE!I  WR, 

CONT  INUE 


1,C0*N 

.EQ.O)  GO  TO  1011 

.GE.100.0.0R.:K!I ,2).GE.100.0)G0  TO  1005 
.GE.100.0.0R.CK!1,4).GE.100.0)GO  TO  1005 
948)P!I.1),P!I,2)  ,P!I.3) ,P!I.4) 

960)CK! I.1),CK(I.2).CK!I t3),CK!I,4) 

948) P! I ,1),PII ,2)  ,P! 1,3) ,P(I,4) 

942)CK!  i,l),CKn,2),CK!I,3)  ,CK(I,4) 

.GE.100.0.aR.CK!l,2).GE.100.0)  30  TO  1015 
.GE. 103.0)  GO  TO  1015 
945)P! I ,1),P!I ,2)  ,P!I,3) 

949) CK!  I.ll.CKll, 2),CK(I ,3) 

945)P!I,1),P!I,2)  ,P!I,3) 

983)CK! i,l),CKlI,2),CK(It3) 


FORMAT  STATEMEMTS 

F0RMAT!1X,A79) 

F0RMAT16F8.5) 

F0RMATI4F8.4) 

F0RMAT12F8.0) 

fOPMAT(9I4) 

F0RMAT!9I3) 

F0RNAT!9I4) 

F0RMAT!F8.3,F8.5, I8,2F8. 1) 
F0RMAT!9I8) 

F0RMATIF8.6,F8.5,  F8.6) 
F0RMATI9I4) 

F0RMAT!F8.6,2F8.5,F8.6) 

F0RMATI3F8.4) 


CLOSE!IWR,IOSTAT=ICK) 

IFIICK  .NE.  0)  WRITE!IQT,921)EEX 
FORMAT! *TROUBLE  CLOSING  PRINTER  OUTPUT 


F ILE*  ) 


lo  OO  OIJ  (-)(_>  Oi 


AeeEMQi^  fi 

Q1  COMPUTER  GENERATED  MODEL  VERIFIER  PROGRAM 


ccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccc 

PROGRAM  Q1 

PROGRAM  Q1  IS  DESIGNED  TO  ACCEPT  VARIABLE 
OIMENSlSNSf  INPUT  A  DATA  FILE  OF  NODE 
TEMPERATURES  GENERATED  BY  TVSSI  AND  CREATE 
AN  OUPur  FILE  OF  HEAT  FLOW  VALUES  ACROSS 
1)  VAPOR  CONDENSATE  INTERFACE.  2)  FIN  AND 
TUBE  OUTER  PERIMETER.  AND  TUBE  INNER  WALL 


C 

C 

DIMENSION  CK(300.4).T(300) 

INTEGER  PI300i4).C0.N.X.ZW.ZWl.ZW2.ZW3.Vl.V2 .W.U. 
double  PRECISION  CTHM .CTHL. HI .HS .RO. RI . A .6. C . D. 
♦VAL4.VAL14.VAL6.VALI6.VAL8.VAL13.VAL5.VAL7. VALID. 
*A2.A3.A4.A5.WINC1 .WINC2.HINC1.HINC2.H1NC3.ALINE. 
♦TIN.TS.T.Xl tX2.VAL2.VAL12.VAL20«Al«RADH.RADL«EL 
CHARACTER  *79  TITLE 
CHARACTER  *79  FNAME 
C 

DATA  IOT.IN.IPR.IWR/6.5.4.8/ 


OPEN(IPR.FILE=»PRN».STATUS=*NEW* .FORM=*FORMATTED' , 
*IOSTAT=ICKI 

IFdCK  .NE.  0)  WRITE!  I0T.20  ) 


FORMAT! 'TROUBLE  OPENING  PRINTER  OUTPUT  FILE* I 
WRITEIIOT. 17) 

FORMAT!/'  ENTER  THE  TITLE  OF  THIS  STUDY.  79  ', 

♦  '  COLUMNS  MAX.  THIS  TITLE  WILL  APPEAR'./,'  AT' 

*  '  THE  TUP  OF  EVERY  PRINTED  PAGE  OF  OUTPUT:') 
REAO!IN.18)  TITLE 

FORMAT!  A79) 

WRITEIIOT,  1) 

FORMAT!/'  ENTER  INNER  RADIUS  OF  TUBE:') 

REAO!IN,2)  RI 

F0RMATIBN.F10.5) 

WRITEIIOT, 3) 

FORMAT!/'  ENTER  OUTER  RADIUS  OF  TUBE:') 
REAOI1N.2)  RO 


WRITEIIOT, 

FORMAT!/' 

*  'FIN  TO  L 
READ!  IN. 2) 


I0T,4) 

!/•  ENTER  DISTANCE  FROM  OUTER 
TO  LIQUI D  SURFACE:  ' ) 


EDGE 


WRITEIIOT, 5) 

FORMAT!/'  ENTER  HEIGHT  OF  FIN') 
REAO!IN,2)  B 

WRITEIIOT, 6) 

FORMAT!/'  ENTER  DISTANCE  BETWEEN  FINS') 
REA0(IN,2)  C 


lOOOOUUOU' 


10 

0 

11 

9 

12 

0 

13 

14 

15 

16 

99 


C 

C 


25 

C 


WRITE(I0T,7) 

FORMAK/*  ENTER  WIDTH  OF  FIN*) 

REA0(IN«2)  0 

WRITEdOTtS) 

FORMATi/'  ENTER  HEAT  TRANSFER  COEFFICIENT  OF 
*•  TUBE  INNER  SURFACE:*) 

REA0<IN«2)  HI 

WRITE(I0Ty9) 

FORMAT!/'  ENTER  HEAT  TRANSFER  COEFFICIENT  OF  *, 

♦  'OUTER  LIQUID  SURFACE:*) 

REA0(INf2)  HS 

WRITEIIOTylO) 

FORMAT!/*  ENTER  THERMAL  CONDUCTIVITY  OF  FINS  AND 
**  TUBE:*) 

READ(IN«2)  CTHM 

WRITECIOT,  11) 

FORMAT! /•  ENTER  THERMAL  CONDUCTIVITY  OF  FLUID  *  , 

*  •  SURROUNDING  TUBE  AND  FINS:*) 

REA0!INt2)  CTHL 


WRITE!I0T,12) 

FORMAT!/'  ENTER  TEMPERATURE  OF  FLUID 
REA0(IN,2)  TIN 


INSIDE  TUBE:*) 


WRITE(I0T#13) 

FORMAT!/'  ENTER  TEMPERATURE  EXTERNAL  TO  *t 
*'  SURROUNDING  FLUID:') 

READ(IN.2)  TS 
WRITE!10T*14) 

FORMAT!/*  ENTER  THE  NUMBER  OF  COLUMNS  DESIRED  ', 

*•  FOR  YOUR  model:  *) 

REA01IN,15)  N 
FORMAT!BNy  12) 

WRITE<I0T»16) 

FORMAT!/'  ENTER  THE  NUMBER  OF  ROWS  DESIRED  FOR  ', 

♦  '  THE  TUBE:*) 

READ! IN. 15)U 
WRITE!I0T,99) 

FORMAT!/*  ENTER  THE  NUMBER  OF  COLUMNS  DESIRED  FOR  *■ 
♦'  THE  FIN:  *) 

REA0IIN«15)ZW 


CO=INT!  300-0/FL0ATIN)  ) 
INITIALIZE  THE  P  AND  CK  ARRAYS 

DO  25  I*1,C0*N 
DO  25  J-1,4 
P!  I. J)«0 
CK!i  fJ)-0.0 
CONTINUE 

VALUES  COMMON  TO  ALL  RUNS 

L21=C0*N 

L22=2 

L23=0 


77 


c 

c 

47 

C 

C 

48 


C 

C 

49 

50 


C 

c 

51 

C 

c 

55 

C 

c 

c 

c 

c 


c 

c 

c 


c 

LOO 

c 

105 

c 

110 

c 


P(I»1J=(I-1I*10+1 
GO  TO  46 


P(Itll=(I-N)*10+l 
P(If2)=(I-l)*10+l 
P(  I#3)=(I-H  )*10+1 
P( It4)=3024 
GO  TO  55 


P(lt2J=(I-li*10+l 
P  (  It  1)=(  I-N)*10+l 
P(  1,31=3024 
GO  TO  55 


P( It2)=( I-l)*10+l 
GO  TO  50 


CONTINUE 

CALCULATE  HEIGHT  AND  WIDTH  INCREME'JTS 

ZW=INT(FL0AT(N)/3«0) 

ZW1=N-ZW-1 

ZW2= INTI (FLOAT (CO  I -10. 01/2.0) 

ZW2=2 

ZW3*C0-(U+1J-ZW2 

WINC1=0/FL0AT(ZW) 

WINC2=C/FL0AT( ZWl i 
H INC 1=(R0-RI)/ FLOAT! U) 
HINC2=B/FL0AT(ZW3) 

HINC3=A/FL0AT(ZH2) 

EL=(  WINC1+WINC2I/2.0 

COMPUTE  THE  RADIUS  OF  EACH  ROW  OF  VJODES 

00  1000  1=1, CO 
IF(I.GT.l)  GO  TO  100 
RADH=RI+0.5*HINC1 
RAOL=RI 
GO  TO  150 

IF(I.GE.(U+1) )  GO  TO  105 
RAOH=RI+(FLOATCI)-0.5)*HINC 1 
RAOL=RI+(FLOAT(II-1.5)*HINC1 
GO  TO  150 

IF(I.NE.(U«-1) )  GO  TO  110 
RADH=R0-»-HINC2/2.0 
RADL=RI+(FL0AT(I)-1.5I*HINC1 
GO  TO  150  ^ 

IF(I.GE.(11+ZW3)I  GO  TO  115 
RAOH=RO+(FLOAT(I l-(FL0AT(U)+0.5) )*HINC2 
RADL=RO+(FLOAT  (n-(FLQAT(U)+1.5)  )*HIMC2 
GO  TO  150 


IF(I  .NE.(U«14-ZW3)  )  GO  TO  120 
RADH=R0+B+Hl'(C3/2.0 
RADL=R0+B-HI'(C2/2 .0 


^  ^  -  V  «^-  -.  . 


GO  TO  150 
C 

120  IFd.EQ.COl  GO  TO  125 

RADH=R0+B+(FL0ATCI )-(FLaAT(U}+0.5)-FLOAr(ZW3) )*HINC3 
RA0L=R0+B+(FL0AT(  IJ-(FL3AT(  U)+l.  5 J -F LOAT  (  ZW3  )  )<'HINC3 
GO  TO  150 
C 

12  5  RADH=R0+B+A 

RA0L=R0+B+A-HI  NC3  /2 .0 
C 
C 
C 

150  ALINE=RA0H**2-RADL**2 
C 

C  COMPUTE  INTER-NODE  CONDUCTANCE  VALUES 

C 

C 

185  IFd.GT.l)  GO  TO  187 
HTC=HI 
R0=RI 
GO  TO  189 

187  IFd.NE.CO)  GO  TO  190 
HTC=HS 
R0=R0+A+B 

189  Al=HTC*WlNCl*R0/2.0 
A2=A1*2.C 
A3=A2*EL/taINCl 
A4=HTC*WINC2*R0 
A5=A4/2.0 

C 

C 

190  IFd  .GE.{U+2)  I  GO  TO  200 
C 

VAL4=RADL*WINC1*CTHM/HINCI 
GO  TO  300 

200  IFd.GE.dU+2)+ZW3J)  GO  TO  250 
VAL4=RA0L*i^  INC  1*C T HM/ HINC2 
GO  TO  300 

250  VAL4=CTHL*RADL*WINC1/HINC3 

300  VAL2=VAL4/2,0 

C 

C 

C 

C 

IFd.GE.(U+lii  GO  TO  310 
VAL14=CTHM*RA0H*W INCl/HINCi 
GO  TO  400 

310  IFCI  .LE.(U+ZW3  )}  GO  TO  3  50 
VAL14=CTHL*RADH*W INC1/HINC3 
GO  TO  400 

350  VAL14=CTHM*RADH*WINC1/HINC2 
400  VAL12=VAL14/2 .0 

C 
C 

c 

IFd.GE.(U+2n  GO  TO  410 
VAL6=CTHM*RA0L*EL/HINC1 
GO  TO  500 

410  IFd.GE.I  (U+2)+ZW3))  GO  TO  450 

VAL6  =  dCTHM*WINCl+CTHL*WlNC2»*RADLI  /  (2.0«HnC2i 
GO  TO  500 

450  VAL6=CTHL*RA0L*EL/HINC3 
C 
C 
C 


80 


500 
61  0 


620 

C 

c 

c 

650 


670 


680 

700 

C 

C 

c 


76  0 


770 

C 

C 

C 

800 


810 


82  0 
C 
C 
C 

845 


850 

C 

C 

86  0 


IF(1  .GE.(U+i) )  GO  TO  610 
VAL16=CTHM*EL*RA0H/HINC1 
GO  TO  650 

IF(I.LE.IU+ZW3 ))  GO  TO  620 
VAL16=CTHL*EL*RA0H/HINC3 
GO  TO  650 

VAL16=(  (CTHM*WINCl+CTHL*WINC2)*RADH  )/  <2.0*HINC2) 


IF(I .GE. (U+2) i  GO  TO  670 
VAL8=RA0L*WINC2*CTH«/HINC1 
GO  TO  700 

IFn.GE.((U+2)+ZW3))  GO  TO  680 
VAL8=CTHL*RA0L*WINC2/HIMC2 
GO  TO  700 

VAL8=CTHL*RA0L*WINC2/HINC3 

VAL10=VAL8/2.0 


IFn.GE.(U+ll  )  GO  TO  730 
V  ALl 8=C  THM*RADH4  M INC2/H1 NCI 
GO  TO  750 

IF(I.LE.(U+ZW3))  GO  TO  740 
VAL18=CTHL*RA0H*WINC2/HINC3 
GO  TO  750 

VAL18=CTHL*RADH*WINC2/HINC2 

VAL20=VALl8/2.0 


IF(I.NE.<<U+IJ+ZW3)J  GO  TO  760 
VAL5*(CTHL*(RA0H**2-(R04BJ**2)+CrHM*( (RO+Bi 
♦♦♦2-RA0L**2))/ (2.0*WINC1 ) 

GO  TO  800 


IF(I  .GT.(<U+1)+ZW3I)  GO  TO  770 
VAL5=CTHM*ALINE/(2.0*WINC1) 


GO  TO  800 


VAL5=CTHL*ALINE/(2.0*WINC11 


IF(I .GE.(U+1) )  GO  TO  810 
VAL7=CTHM*ALIN£/C  2.0*WINC2) 

GO  TO  845 

IF<I.NE.(U+1)  »  GO  TO  820 

VAL7=(CTHL*(RADH<‘*2-R0**2  }+CTHM*(R0**2-R A0L**2)  ) 
«/ (2.0*WiNC2) 

GO  TO  845 

VAL7=CTHL*ALlNE/( 2.0+WINC2) 

COUPLE  CONDUCTANCE  VALUES  TO  CONNECTION  TERMS 


IF(I.EQ.l)  GO  TO  905 
IF(I.EQ.CO)  GO  TO  912 
Q=N*n-l) 

DO  900  J=Q+lfQ+N 


IF(J.EQ.(Q+1I  )  GO  TO  875 
IF(J.LT.(a+ZH  +  l))  GO  TO  880 
IF(J.EQ.(Q*Zrt*H)  GO  TO  885 
IFIJ.NE.(Q+Nn  GO  TO  893 


C 


c 


c 

c 

880 


C 

C 

885 


C 

C 

890 


C 

C 

875 


C 

900 

901 
C 

c 

905 


C 

C 

906 


C 

C 

90  7 


C 

c 

908 


CK(J«1}=VAL13 
CK<J,2»=VAL7 
CK(J,3)=VAL20 
GO  TO  900 


CK( J,1)=VAL4 
CK(J,2i=VAL5 
CK(  J,3I=VAL5 
CK(J,4)=VAL14 
GO  TO  900 


CK(Jtl)=VAL6 
CKi J,2)=VAL5 
CKfJ,3)=VAL7 
CK(J,4J=VAL16 
GO  TO  900 


CK(Jtl)=VAL8 
CK< J,2)=VAL7 
CK(J,3)=VAL7 
CK<  J#4i=VAL18 
GO  TO  900 


CK(J,1)=VAL2 
CK(J,2)=VAL5 
CK(  J,3J=VAL12 

CONTINUE 
GO  TO  1000 


DO  910  J=1,N 

IF(J.GT.l)  GO  TO  906 

CK(Jtli=VAL5 

CK(J«2J=VAL12 

CK(J,3I=A1 

GO  TO  910 


IFU.GE.ZW+11  GO  TO  907 

CK(J, 1I=VAL5 

CK( J,21=VAL5 

CK(J,3J=VAL14 

CKI J,4J=A2 

GO  TO  910 


IF(J.NE.(ZM-Mn  GO  TO  908 

CK(J,1J=VAL5 

CK(J,2)=VAL7 

CK(J«3)=VAL16 

CKIJ,4)=A3 

GO  TO  910 


IFIJ.EQ.N)  GO  TO  909 

CKU,1)=VAL7 

CKI  J,2)  =  VAL7 

CKIJ,3J=VAL18 

CKI  J,4»=A4 


82 


_ « 


C  READ  IN  NODE  TEMPERATURES  GENERATED  BY  TVSSI 

C 

C 

DO  924  J=l»25 

READ! 16,1045) (TII ) , I=< ( J-1) ♦12+1 ) , J*12 ) 

92  4  CONTINUE 


C 

92  5 
C 

926 


92  7 


92  8 
C 
C 
C 


c 

c 

949 


950 

955 


96  0 


965 

9000 

970 

975 


98  0 


2000 


C 


MRITE(I0r.925) 

FORMAT!/'  WRITING  FILE  TVSSI') 

M  R IT  E ( I OT , 926 ) CO^N .CO , N , ZW2 
FORMAT!/'  THIS  ',13.'  NODE  MODEL 


ruivnMi,/'  inid  nuuc  nuuci.  i\u 

♦'  AND*  ,12,'  COLUMNS.  THERE  ARE  ',12,*  ROWS') 
WRITE!I0T,927)ZW3 

FORMAT!'  IN  THE  CONDESATE  1.AYER  AND  ',12, 

IN  THE  FIN.  THE  FIN  CONTAINS') 
WRITE!IOT,928)ZW 
FORMAT!I2,'  COLUMNS.  ') 

COMPUTE  HEAT  FLOW  ACROSS  DESIGNATED  SURFACES 

QFL=0.0 

QFS=0.0 

QFU=0.0 

QU=0.0 

U0=0.0 


ROWS 


00  2000  1=1, N 

IF!I.GT.!N-ZW-1))  GO  TO  955 

V1=(U+1)+N+1+ZW+I 

V2=( U4N)+1+ZW+I 

IF!V2.EQ.  (U+1)4N)  GO  TO  950 

QFL=UFL  +  (T! VI )-T!  V2) ) ♦CK ! V2 ,4) 

GO  TO  955 

QFL=UFL+1T!V1 )-T ! V2) ) ♦CK ! V2 ,3) 

IF(I.GT.ZW3)  GO  TO  960 
W=N+(U+I )+2+ZW 

QFS=QFS+!T!W)-T(W-1))^CK!W, 3) 

IF!I .GT.IZW+1) )  GO  TO  970 
Xl=!U+ZW3)^N+I 
X2=! U+ZW3+1)^N+I 
IFlI.EQ.l)  G3  TO  965 
gFL=QFU+!TlX2)-T! X1))^CK!X1 ,4) 

WRITE!6,9000iQFU, T !X2 ) ,7 ! XI ) , CK! XI, 4 ) , XI , X2 
GO  TO  970 

QFU=QFU+(T!X2 )-T ( X 1 ) ) ♦CK 1 XI , 3 ) 

WRITEi6,9000)OFU, T!X2) ,T!X1  ),CK!XI, 3) ,X1,X2 
FORMAT!1X,4F10.5,2I3) 

IF!I.GT.1)G0  TO  980 
Q0=00+!T!I)-TIN)4CK!I,3) 

QU=QU+!TS-T 1I+!N^!C0-1)) ) )^CK ! ! I +! N^ ( CO- 1 ) ) ) ,3) 
GO  TO  2000 

IF!I.EQ.N)  GO  TO  975 
QO=QO+!T! I)-TIN)^CKII,4) 

QU=QU+!  TS-T!  I  +  !N^  (CO-D)  )  ) ♦CK  ( ( I +  !  N^! CO-1 ) )  )  ,4) 
CONTINUE 
CIT1=QFU+QFS 
QT0=QT1+QFL 
WRITE  STATEMENTS 
WRITE!IWR,1030) 

WRITE!IWR,1031)RI , RO, B 
WRITE1IWR,1032  ) 

WRITEIIWR.1033)  C,0,A 


c 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1040 

1041 
10  42 

1043 

1044 

1045 

1047 

1048 


WR1TE(IWR»1034) 

WRlTEnWR,1035)TlNtTS 

MRITE(IWRtl036) 

WRITE(lWRtl037)CTHM»CTHL»HI  »HS 
MRITE(IWRtl047iQFU 
WRlTE(IURtl048}QFS 
WRITEUMRf  1040)QT  1 
'rtRiTE(lWRtl041}QF  L 
WRITE(IWR«1042)QT0 
MRITE(IURtl043)QU 
WRITE(IWR.1044)QQ 


FORMATi/*  PIPE  INNER  RADIUS  PIP 

*•  FIN  HEIGHT  M 
F0RMAT(4X,F8.5f'  M  •  ,1 IX,  F  8.  5,  •  M', 
FORMAT!/'  FIN  SPACING  FIN 

♦*  CONDENSATE  HEIGHT  •) 
F0RMAT(1X,F8.5,'  M * ,10X, F8. 5, *  M*. 
FORMAT!/'  TEMPERATURE  INSIDE  PIPE 
*•  OUTSIDE  CONDENSATE  LAVER  •) 
F0RMAT17X,F4.1,'  C  '  ,34X,  F4.  1, '  C) 
FORMAT!/'  K  METAL  K  LIQUID 

♦'  H  IN  CONDENSATE  LAYER') 


PIPE  OUTER  RADIUS 


M' ,5X,F10.5, '  M'  ) 
FIN  WIDTH 

M',7X,F9.8,'  M') 

PE  TEMPERATURE 


FORMAT!/'  K  METAL  K  LIQUID  H  INSIDE  PIPE 

♦'  H  IN  CONDENSATE  LAYER') 

F0RMATIF5.0,'  W/M  •  ,4X  , F5 . 1,  •  W/M' ,3X,F  10.0, '  W/M**2'  , 
♦  5  X,  10.0,'  W/M**2') 

FORMAT!/'  Q  IN  ON  FIN  UPPER  AND  SIDE  SURFACES  ', 


FORMAT!/'  Q  IN  ON  FIN  UPPER  AND  SIDE  SURFACES  ' 
♦F9.3,'  WATTS  ') 

FORMAT!/'  Q  IN  ON  PIPE  SURFACE  ',15X,F9.3,'  WATTS') 
FORMAT!/'  TOTAL  Q  IN  THROUGH  PIPE  AND  FIN  SURFACE  ', 

♦  F5.3,  '  WATTS  '  ) 

FORMAT!/'  Q  IN  THROUGH  CONDENSATE  ',12X,F9.3,'  WATTS') 
FORMAT!/'  Q  OUT  THROUGH  PIPE  INNER  SURFACE  ',3X,F9.3, 
*'  WATTS  ') 

FORMAT!  12F6.2  ) 

FORMAT!/'  Q  IN  THROUGH  UPPER  FIN  SURF ACE • , 9X , F9 . 3 , 

♦  '  WATTS  •) 

FORMAT!/'  Q  IN  THROUGH  FIN  SIDE  SURFACE  •,9X,F9.3, 

♦  '  WATTS  •) 


1051 

C 

C 


CLOSE!IWR,IOSTAT= ICK) 

IF!ICK  .NE.  0)  WRITE!  I0T,1051)EEX 
FCRMAT!'TROUBLE  CLOSING  PRINTER  OUTPUT 


STOP 

END 


FILE*  ) 
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